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We present a unified account for tlie coupled single-hole- 
and spin-dynamics in the spin-gap phase of dimerized and 
frustrated spin-chains and two-leg spin ladders. Based on the 
strong dimer-limit of a one-dimensional ti,2,3- Ji,2,3-model a 
diagrammatic approach is presented which employs a map- 
ping of the spin-Hamiltonian onto a pseudo-fermion bond- 
boson model. Results for the single-hole spectrum are de- 
tailed. A finite quasi-particle weight is observed and studied 
for a variety of system parameters. A comparison with exist- 
ing exact diagonalization data is performed and good agree- 
ment is found. 
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I. INTRODUCTION 

Unconventional quantum-magnetism in low-dimensi- 
onal transition-metal compounds has received consider- 
able interest recently due to the discovery of numerous 
novel materials with spin-i moments arranged in chain, 
ladder, and depleted planar structures. Among these 
compounds are antiferromagnetic chain-systems which 
are intrupsically dimerized, in particular, (VO)2P207Ej, 
CuWoJ, and Cu(N03)2 • 2.5H20a Moreover, quasi- 
one-dimensional (ID) materials have been discovered 
which display^a temperature dependent dimerization, 
eg. CuGeOsQ and a'-NaV2 05El where the former is 
the first igarganic spin-Peierls chainEI and magnetically 
frustratedQiS while the latter is a ^-filled two-leg ladderB. 
Spin-lacLders, both inorganic, eg. SrCu2 03lllj aud, 
CaV205M as well as organic, eg. Cu2(C2Hi2N2)2CUp^ 
have been investigated. Quite recently SrCu2(B03)2li3'Ej 
has beea shown to realize the two-dimensional (2EL) 
versionEjof the ID frustrated Majumdar- Ghosh modelEj. 

In contrast to conventijOnal ID spin-chain materials, eg. 
Sr2Cu03 and SrCu02EZfEa, which display algebraic, al- 
most long-range, spin-correlations and gap-less magnetic 
excitations the new materials are spin-liquids (or dimer- 
solids) with short-range singlet correlations and a gap in 
the spin spectrum. The spjijB^ap phenomenon, hat>Jarcn 
attributed to dimerizationollj and frustrationlijElJ'Eil in 
chain-systems, and it can be interpreted accordingly in 
ladders due to their tpapiogical equivalence to frustrated 
and dimerized chainsE3E3. 

Apart from magnetic excitations the dynamics of 
electronic carriers doped into quasi- ID spin-liquids is 
of interest, in particular because of the discovery of 
superconductivityriL in the two-leg ladder-compound 
Sri4_2:Ca2:Cu2404iE3. Important progress has been 



achieved regarding the uairing correlations and the phase 
diagram of spin-laddersa. Howeyp, only a restricted set 
of primarily numerical studie£a~E3 has focused on the 
spectral properties of single-hole excitations in quasi- ID 
spin-liquids at low doping. In this regime, and due to 
the spin-gap, a quasi-particfe pictureLj, rather than a 
Luttinger-liquid despdptiontJ, is believed to be relevant 
at low energy scalesEII. 

In this paper we detail a theory of single-hole excita- 
tions at half-filling for a dimerized and frustrated spin- 
chain. Particular emphasis will be on the coupling be- 
tween spin and charge degrees of freedom. The paper is 
organized as follows: first we introduce the ii,2,3->/i,2,3- 
model for the dimerized and frustrated spin system. Sec- 
ond we summarize a bond-boson-method to evaluate the 
spin-excitations at half filling starting from the strong 
dimer limit. Third, we map the ti,2,3-t/i,2,3-iiiodel in the 
single-hole sector onto a coupled boson-fermion-model. 
Next, the single-hole excitations of the latter model are 
evaluated by a selfconsistent diagrammatic technique. 
Finally we present results for the spectral properties of 
the single hole, both for ladders and chains, and com- 
pare our findings to existing numerical analysis of finite 
systems. 

II. THE fi,2,3-Ji,2.3-MODEL 

For the remainder of this work we focus on systems 
which are of a Mott- or charge-transfer type of insula- 
tor at half filling and allow for an approximate mapping 
to a i- J-model. To describe the case of either dimeriza- 
tion and frustration of the magnetic exchange coupling 
along the chain or the topology of a two-leg ladder we in- 
voke a i- J- model which includes hopping- and exchange- 
integrals up to third-nearest-neighbors (NNN), i.e. ii,2,3 
and Ji,2.3, sls shown in figure |l|. 
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FIG. 1. ti,2,3- Ji,2,3-model. n labels the dimer bonds. 1 
and 2 refer to non-equivalent sites per unit cell. 



1 



Depending on the case of a ladder or a chain the real- 
space structure of the system is either identical to that of 
figure with J2 = ^2 = or can be obtained by deforming 
the figure such as to place sites 2,n half-way in-between 
sites 1 , n and l,n + I. The Hamiltonian reads 



H = Ht + Hj 



Hj — Ji Si „S2_„ + J2 S2,„Si^„_|_i 

n n 
+ -^3 ^^(Sl,„Si,„+i + S2,.,iS2,„+l), 



(1) 



(2) 



n.{7 n,a 
n,cr 

where Si_„ are spin- 1/2 operators on site i of dimer- 



bond n and 



a(t) 
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,n,~a)y^' are projected 
fermion-operators of spin a on site i, n. Dimerization 
and frustration of the spin-system is expressed via the 
parameters 5 and a where Ji(2) = ■AjII + {~)^) and 
a — J^/Jq. Similarly the hopping integrals are related 
through ti(2) = to{l + (— )(5) and a = ts/to- Ladders are 
characterized hy S = 6 = 1. Regarding the spin-part, 
i.e. Hj, the ground state is known to be a product-state 
of singlets located||-©p the dimer-bonds for a and S on 
the disorder-lineE3ll3 S = 1 — 2a, moreover, the spin- 
spectrum is j gap- .less only at 6 = and < a < ac with 
ac — 0.2411c3c]l. For magnetic couplings which arc me- 
diated by super-exchange and for systems with only a 
single relevant on-site Coulomb-energy scale we expect 
that (1 -I- (5)^ (H- 5) and a for \S, a\ <C 1. 



III. UNDOPED SPIN-SYSTEM 

The spin-part, i.e. Hj, of Hamiltonian (|^) allows for 
a mapping onto a model of hard-core 'bond'-bosons orig- 
inally employed in the context j-of the 2D Heisenberg- 
model for the high-T^ cupratea23. Here we briefly re- 
state the essential features of this mapping. Consider any 
two spin-1/2 operators Si,,i and S2^„. The eigenstates 
of the related total spin are a singlet and three triplets 
These can be created out of a vacuum |0) by applying 
the bosonic operators sjj and tj^ „ with a = x,y, z 



4lo) = ^(in>-UT». 
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4,J0) = -^(|TT>-|U». 
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(ITT> + IU». 
(in> + UT». 



(4) 



where the first (second) entry in the kets refers to site 1(2) 
on dimer n of fig. |l|. On each site we have [s,s'l'] = 1, 
[s^l'^i^t)] ^ 0, and [t„ ^t^] = S^p- The action of Si,„ 
and S2,n in this space leads to the representation 



S? 



(5) 



Here ea/j-y is the Levi-Civita symbol and a summa- 
tion over repeated indices is implied hereafter. The 
upper(lower) subscript on the Ihs. of (^ refer to up- 
per(lower) sign on the rhs.. The bosonic Hilbert space 
has to be restricted to the physical Hilbert space, i.e. to 
either one singlet or one triplet per site, by the constraint 
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(6) 



The representation (^) can be inserted into (|^) yield- 
ing an interacting-tjose-gas Hamiltonian accompanied by 
the constraint (|6|)E3. This Hamiltonian is diagonal intra- 
dimer-wise and contains two-particle interactions which 
are only of inter-dimer type. At the point of complete 
dimerization, i.e., i5 = 1 and a = 0, the inter-dimer 
interactions vanish leaving a sum of purely local dimer 
Hamiltonians, each of which has a singlet ground state. 
This renders the global ground-state a product of sin- 
glets localized on the dimers with the excitations being 
a set of 3^-fold degenerate triplets. Off the dimer-point 
the inter-dimer interactions can be treated approximately 
by a linearized Holstein-Primakoff (LHP) apptoach, de- 
tails of which can be found in the literaturcEZrc^l. The 
LHP method retains spin-rotational invariance and re- 
duces Hj to a set of three degenerate massive magnons 
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The '[] '-bracketed sign on the rhs. in ( pi] ) refers to the 
quantity v on the Ihs.. The spin-gap A = minjwfc} re- 
sides at k = 7r(0) with A = ^Jf - (+)Ji(2J3 - J2) for 
2J3 > (<) J2- Note that because of (||) the ground state 
\D), which is defined by ja,k\D) = 0, contains quantum- 
fluctuations beyond a pure singlet product-state. 

To leading order the dispersion Wfc is , identical 
to perttujbative expansions, both, for chainsoc3 and 
Iaddergl3e3 and it has been applied to model ip^. 
elastic neutron j-scattering (INS) data for CuGeOacHi 
and a'-NaV205E3. Beyond the LHP approach triplet- 
interactions lead to a renormalization of Uk and to the 
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formation of multi-magnon bound statesc3c3. More- 
over the constraint although a hardcore repulsion, 
has been considered perturbatively .yelyi pg on the case 
of low triplet-density at A/ J < lliaSEl. These renor- 
malizations of the triplet dispersion will be discarded in 
our evaluation of the single-hole spectra, in particular 
since their dominant effects can be accounted for semi- 
phenomenologically by adjusting the size of the LHP 
spin-gap. This is of significance when compari ng si ngle- 



hole spectra to numerical studies, as in section VII 



IV. SINGLE-HOLE HAMILTONIAN 

In this section we map the ti_2,3-</i,2,3-IIamiltonian of 
(|l[^) in the single-hole sector onto a model of pseudo- 
fermions interacting with the bond-bosons. To this end 
we note that a dimer-bond occupied by a single hole can 
be labeled by introducing an additional pseudo-fermion 
(holon). I.e., instead of the bond being in one of the 
states given by it can also be in the state 



(12) 



where the l.h.s. denotes the vacuum, i.e. |0), with a 
single dimer-bond at site n in a one-hole state of spin a 
with j = 1,2 referring to the two positions available to 
the hole on the bond. The operators ^ ^ are required to 
obey fermionic anticommutation relations. In the single- 
hole sector each dimer-bond can only be in exactly one 
of the states given by (Q) or (|l2|) . Therefore an extended 
hard-core constraint has to be satisfied 



i=i,2 



(13) 



with a summation over repeated spin-indices implied. 
Creation of a single physical hole in the (half-filled) 
ground state \D) of the spin-system is achieved by ap- 
plying the two-particle operator 



V2 J'"^'^ 



(14) 



where Pj ^ +(-), j = 2(1) for j = 1(2) and pa = +(-), 
a =1 (t) for a —"f (J,). In terms of the previous equation 
the creation of a physical hole can be interpreted has 
the removal of a spin-dimer followed by the creation of 
dimer hole-state, where the particular linear combination 
of the „ g., s„, and t^, „ operators on the r.h.s. ensures 
that the total spin is 5 = 1/2 and = ±1/2. Using 
the constraint (|l3|) it is straightforward to show, that the 
r.h.s. of (|^ indeed satisfies the usual Hubbard-operator 
algebra. 

Inserting (|lj) into we have to distinguish be- 

tween various cases, i.e., (i) intra-dimer hopping, (ii) 



inter-dimcr hopping, and (iii) exchange-scattering. Inter- 
dimer hopping can be either spin-diagonal or accompa- 
nied by spin-flip scattering. This includes (ii,a) singlet- 
singlet, (ii,b) singlet-triplet, and triplet-triplet transi- 
tions of the spin-background upon hole-hopping, the lat- 
ter may occur with a change in the spin-quantum num- 
ber of the background of either (ii,c): ASz = or (ii,d): 
AS, = 1. 

Processes of the type (i) and (ii,a-d) result from direct 
insertion of ( p^ into (^). To express the inter-dimer 
exchange-scattering, i.e. process (iii), one has to realize 
that the spin operator on a dimer in a one-hole state 
needs to be expressed in terms of the pseudo-fermions 
rather than the bond-bosons, i.e. using 



(15) 



and setting j ~ i and m — n, where To-jg-g ^-re the Pauli 
matrices and (SJ^^)^ ~ ^nm- Substituting this repre- 
sentation for either a right or a left spin-operator of the 
inter-dimer part of the exchange into (H) leads to single- 
hole exchange-scattering terms cx Ji,2.3- Summing all 
contributions we arrive at the Hamiltonian 

H = -hY^ al,n,aa2,„,^ + h.c. 

+ Y 4,n,^ai,n-l.<T + ^-C- 

+ hsY^i (s2Vi,„ - si\^) + h.c. 

n 

+ ~2~'^^^'^ (^n-l,n-l ~ ^n+l.n+l) + ^-C 
n 

+ t,s Y tU-iy-' (s^^i,„ + sf^i,„) + h.c. 

j,n 

2 E] '^n~l'^n'A,n,a-'^l,n-l,(7 + 
n,(T 



J2 

2 
2 



9 5Z * ^ + ^n+l,n+l) 

n 

^ (t,\ X t„) (S^ii,„+i + ^lUn-l) ■ (16) 
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with tjj = itx.m'ty,m'tz.7i)- This Hamiltonian is spin- 
rotationally invariant which is consistent with the ground 
state displaying no long-range magnetic order. According 
to the LHP approximation we have replaced the singlet 
operators by C-numbers, i.e. s. Within the LHP ap- 
proach s = 1 to lowest order. Regarding the single-hole 
Hamiltonian we improve upon this by assuming that the 
condensate density of the singlet is determined by satis- 
fying the hardcore constraint (|^) on the average, i.e. by 
setting s„ = = (s„) = s with 



(17) 



In principle this equation can be used to determine the 
LHP magnon-dispersion selfconsistently which however 
will not be done here. For the particular case of a lac 
i.e. J2 pi^2 ~ 0, Hamiltonian (|l^) agrees with refs.l 
See ref.EJ for a related representation of the bilayer HT^- 
cuprates. 



V. SINGLE-HOLE GREEN'S FUNCTION 

In this section we evaluate the retarded pseudo- 
fermion, and physical-fermion Green's function 

GAkA) = ~ie{t){D\{^,Jt),4j\D) (18) 
Glik, t) - -tQ{t){D\{cj,^Jt), c^ljlD), (19) 



where VL = {al 



2,k,i 



and 



cr' '-2,k 



i.e. both, Ga{k,t) and G%{k,t) are 2x2 matrices. We 
proceed via standard diagrammatic techniques to to eval- 
uate the Green's functions. To this end approximations 
have to be ma|ie. First, and in the limit of large dimer- 
ization we fincE3, that renormalization effects due to the 
two-triplet vertices in (^6|) are of minor importance and 
we will discard them in the following. Fourier transform- 
ing the corresponding simplified Hamiltonian and intro- 
ducing the Bogoliubov representation of the triplets, i.e. 
(|), we obtain H = Hq + V with 



k,(7 

^ fe,g,CTl,CT2 

where 



Ek = 
with 



l^k 



Mkq = 



2,kq 



(20) 



(21) 



Mfc = ■s^is cos k 



D(q, CO') 
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G(j(k-q,CO- Cij) 



FIG. 2. SCBA self-energy. 



mi^kq = Uqs{—e'"^ - ts cos {k - q) 



■ cos q) 



+Vqs{^ — ts cos k 
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— cosq) 



m2,kq 



t2 
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ik \ 



(22) 



The structure of (pO|-22) is reminiscent of the Hamil- 
tonian which occurs in the pseudo-particle description 
of the single-hole dynamics in an antiferromagnetic 
(AFM) spin-background which has been studipd exten- 
sively in the context of the HTc-compoundsL3~E3. These 
studies have shown the so-called selfconsistent-Born- 
approximation (SCBA) to the pseudo-fermion Green's 
function to be in satisfactory agreement with informa- 
tion available from exact-diagonalization (ED) studies on 
finite lattices. The SCBA amounts to a non-crossing, 
infinite order resummation of fermion-boson scattering- 
diagrams to the pseudo-fermion self-energy, neglecting 
vertex corrections as depicted in fig. 2. In an AFM 
spin-background the neglect of vertex corrections caii Jbe 
justified in the limit of high coordination numberE3~EZl. 
Regarding a dimer-state, i.e. |-D), we are unaware of 
a similar simplification. However, at the one-loop level 
we have checked numerically that the inclusion of-vjcrtex 
corrections leads to relatively minor changes onlyEa. 

The SCBA self-energy Eo-(fc,a;) of figure 2 is evalu- 
ated in terms of the bare pseudo-fermion Green's func- 
tion G"(fc,w) with respect to Hq of (20), i.e. 



G^ik^Lo)^ 



(w - fik)^ - |efcp 



) (23) 



and the retarded triplet Green's- function Da{k,t) — 
-iQit){D\[-f^^^{t),-fl^J\D) the Fourier transform of 
which follows from Eq. (^ 

D^{k,Lu)^^—. (24) 



where it is assumed that uj = ui + irj with 77 — > 0+ in 
(^,^J) as well as in all other retarded propagators in the 
reraainder of this paper. Inserting (^ 24) into fig. 2 we 



gell 



E,(fc, ^) - ^ E (fc -q,u;- Uq)-^ 



E,ik-q,w-Wg)]-'Ml 



kq 



(25) 



from which the pseudo-fermion Green's function is ob- 
tained as usual via G^(fc, w) — [Gj^(fc, W)~^ — Ecr(fc, W)]~^. 
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FIG. 3. Structure of the physical Green's function 
G%{k,u!). Solid, dashed, and dotted lines refer to 
pseudo-fermion, triplet, and singlet Green's functions. 
Dashed box denotes reducible two particle vertex. In LHP 
approximation dotted lines collapse to unity. 
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using the RPA as shown in fig. A comparable anal- 
ysis of the effects of ground-state quantum fluctuations 
on the single-hole spectrum has been carried out in jOn 
AFM spin-background for the square-lattice t- J-modelS. 
Within the RPA the physical Green's function is given 
by 

ak,aG„{k, W)ak,a- + ^ Pkq.aGa{k -q,UJ - UJq)Pkq.a 
1 

( 1 



-G„(k,uj) 



§MkqG,{k-q,u;-u;q)Ml^ 
1 



'k,a 



M 



kq, a Pkq, a Pkq, o 
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kq, a 



"■k,a 



+ 4- 



\ 



Pkq, tJ M|jq ^ M|jq o Pkq, a 

FIG. 4. LHP+RPA approximation to the physical Green's 
function G'a-{k,ui). 



The excitations of physical significance, eg. in a pho- 
toemission experiment, areJJjie creation of holes via the 
two-particle operators (pk.a-^- The corresponding two- 
particle Green's function is depicted in fig. ^ where 
the vertices ak,a and (3kq,a correspond to the different 
ways by which the physical hole can be projected onto a 
pseudo-fermion and a dimer spin-state, i.e. a singlet or 
a triplet. This projection accounts for the ground-state 
quantum fluctuations of the spin system. 



ak,a = s{D\lpk,a(l)k.a\D) = ^ ^ 



(26) 



/3feg,<T = {D\{jq^l:k^q)i^(l)k,a\D) = T]J ^^q'^^ (27) 

where, according to the LHP approximation, the singlet 
creator has been replaced by a C-number and 



J2q MkqGa [k - q,UJ - UJq)(3kq., 
Ga{k,Uj)ak,cr 



(29) 



where the element-by-element multiplication of the ma- 
trix and vector entries in the matrix-product on the last 
three lines of this equation are 2x2 matrices operations. 
To arrive at (|2|) the LHP dynamics of the bond-bosons 
has been used which implies a constant factor of unity 
only in case of the singlet propagator. 

Concluding this section we note the sum rule 



(30) 



for the exact physical Green's function with A'^{k,uj) — 
— ifm G^{k,uj) at half-filling which differs from that for 

free fermions since the c\^^ ^ are Hubbard operators. In- 



tegrating (30) we find that the combined LHP and RPA 
approach is consistent with this sum-rule yielding 



) 

2 [O 



N 



where the term in the last bracket is one due to (1 



(7,Vfc-9)iii = 

■^{Tjz..qiJk-q,l{n + ilx,q T ily ,q)^k-q.] {[)) (28) 

refers to the proper linear combination of spin- 1/2 and 
spin-1 which ensures that the hole created by ipk^a has 
spin S = ^ with Sz — In contrast to the vertices 

( p6| , p7| ) their respective products and squares are spin- 
independent. 

For a finite albeit small density of triplet-pairs in the 
ground state, i.e. < i;, -C 1, we proceed by a pertur- 
bative evaluation of the two-particle propagator of fig. 



VI. DIMERIZED CHAIN SPECTRA 

In the remaining sections of this work we will detail 
various results of the numerical solution of the SCBA and 
RPA equations. We begin with the dimerized-chain limit 
of the ti,2,3-</i,2,3-inodel, i.e. with J3 = is = 0. In view of 
the non-zero dimerization it seems convenient to visualize 
the lattice geometry as that of a linear chain with a unit 
cell containing two electrons and a lattice constant 2a, 
set equal to unity, where a is the inter-site distance. In 
contrast to this, existing ED studies frequently employ a 
different representation in which each unit cell contains 
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a single site only. This representation can be obtained 
by introducing the fermion operators 

d| = -i= (^Cl,k,a + e''^C2,k,a^ , (32) 

which we term d-electrons hereafter. The d-Green's- 
function G'^{k,uj) reads 

G^(fc,c.) - i[(G^(fc,c.))n + (G^(fc,c.))22 

+ e-^^(G^(fc,c.))i2 + e'^(G^(fc,c.))2i], (33) 

with A'^{k,uj) ~ —-Im G'^{k,uj) beeing the correspond- 
ing spectral function. 

Selfconsistent solutions for the SCBA can be obtained 
quite easily by numerical iteration on fairly large lattices 
and for small, but finite imaginary broadening 77 ^ 1. 
Convergence of the iteration is achieved within a few, 
typically 3-20, cycles depending on the size of the spin 
gap. Moreover, finite size scaling analysis can be used 
to determine an approximate system size at which fur- 
ther increase in the number of lattice sites leads to no 
additional change in the pseudo-fermion Green's func- 
tion at fixed 77, i.e. the thermodynamic limit. Inserting 
SCBA Green's-functions from this limit into the RPA 
( p9| ) we obtain typical spectra as displayed in fig. ^. The 
number of sites is 512 = 2 x N where TV is the num- 
ber of dimers. The dimerization of the hopping integrals 
and the exchange couplings have been chosen indepen- 
dently with ti = t2 and 6 such as to result in a spin 
gap of A/ti = 0.2Jo/ii. Since Ji = Jo(l ± S) one finds 

A/h = [2 ((52 + S)]^/^jQ/ti. Figure |a) (b)) refers to 
Jo/ti = 1(0.5) and A/ti = 0.2(0.1). All spectra are 
particle-hole inverted, i.e. uj —^ — w, such as to place the 
first electron-removal state at the lowest binding-energy. 

The figure displays two dominant dispersive regions 
which are related to the first three terms in Hamilto- 
nian (p^). These terms allow for coherent hole- motion 
without triplet excitations of the spin-background and 
lead to two tight-binding-type bands which are accounted 
for by Hq in Hamiltonian (|2^). These tight-binding 
states are renormalized and broadened upon coupling 
to multi triplet-excitations via the SCBA. In particular 
the states of high binding-energy are broadened quite 
strongly. Comparing fig. Ha) vs. b) it is obvious that 
the amount of the spectral redistribution depends on the 
relative size of Jo/^i, where decreasing Jq/^i implies a 
larger incoherent part of the spectrum. Additionally the 
dispersion flattens as Jo/ti decreases. 

Finite-size analysis of the two sharp structures at the 
high-energy edge of the spectrum in figure ^) (b)) re- 
veals that for any finite spin-gap these structures refer 
to poles of the Green's function on the real axis which 
remain separated by a gap from the continuum of multi- 
magnon shake-offs. The spectral weight of these poles re- 
mains finite in the thermodynamic limit. Therefore, the 
first electron removal state is of a quasi-particle type. 
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FIG. 5. Spectral function of the dimerized chain for var- 
ious values of momentum fc = to tt and ti = t2, N — 256, 
and rj = 0.05. (a) Jo/ti = 1, A/ti = 0.2 and (b) Jo/ti = 0.5, 
A/ti = 0.1. 

The first of these two poles, i.e. at higher binding- 
energy, can be traced back to a quasi-particle excitation 
in the pseudo-fermion SCBA-spectrum which is shad- 
owed in the physical fermion Green's function. The sec- 
ond pole at lowest binding energy, i.e. the first elec- 
tron removal state, arises from a zero of the determinant 
of the RPA-matrix propagator in (|2^). This can be in- 
terpreted in terms of a bound state between the afore- 
mentioned pseudo-fermion quasi-particle and the ground- 
state quantum fluctuations of the spin-system. 

Figure ^ refers to the dependence of the spectra on the 
spin-gap displaying two additional values of S different 
from that used in flg. ||. We find that the gap between the 
first electron-removal state and the second pole increases 
upon increasing the spin gap. Moreover the dispersion 
decreases for smaller values of S. For A/ti — 0.05Jo/ii 
and Jo — ti, fig. ||a.), an intense low energy band arises. 
The dispersion of the two dominant spectral regions at 
5 <C 1 are reminiscent of similat-findings for the infinite- 
U Hubbard-chain at half fiUingEj. 

Figure ^ shows the weight, i.e. the Z-factor, of the 
quasi-particle pole at wave vector fc = 7r/2 as a func- 
tion of various parameters, S, Jo/ti, and Js/ti. The Z- 
factor has been determined by fitting a Lorentzian to 
the quasi-particle peak. The renormalization effects are 
quite strong for those parameters we have considered as 
Z is reduced substantially from the non interacting value 
of one half (|30|). This figure demonstrates that a finite 
spin-gap stabilizes the quasi-particle excitations, while 
approaching the Heisenberg point, i.e. 5 = 0, J3 = 0, 
as in fig. ^) suppresses the Z-factor due to the de- 
crease in energy of available triplet excitations. In fig. 
the dimerization is kept at a constant value to result 
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FIG. 6. Spectral function of the dimerized chain for various 
values of momentum fc = to vr and ti = t2 = Jo, N = 256, 
and Tj = 0.05. (a) A/ti = 0.05 and (b) A/ti = 0.7. 
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FIG. 7. Quasiparticle weight Z{k = 7r/2) as a function (a) 
of A/fi for Jo = ti, J3 = 0, (b) of Jo/h for A/ti = 0.2 Jo /h, 
i.e. 5 = const., J3 = 0, and (c) of Js/ti for 5 as in (b), 
Jo =ti. N ^ 256, ti = 

in A/ti = 0.2Jo/ti with only Jo/ii varying. This leads 
to almost no change in the Z-factor. Increasing frustra- 
tion, as in fig. ^) increases both, the average absolute 
value of the vertex- function Mkq of (^ij) and the spin-gap 
A. These effects compete leading to only a slight increase 
of the Z-factor as a function of J3. Here S is still kept 
at the constant value of fig. 0b). While the single-hole 
case does not resemble a finite hole concentration it is 
tempting to relate the finite Z-factor in the presence of a 
spin-gap to a Luther- Emery-liquidcJ scenario in contrast 
to the Luttinger-liquicS at vanishing spin-gap. 

Next we turn to a comparison of our diagrammatic ap- 
proach with results obtained from exact diagonalization 
(ED) of finite chains. To this end, we contrast a single- 
hole spectrum for S — 0.048 and ti — t2 — Je-ircpro- 
duced from the work of Augier and coUaboratorsEil in fig. 
^) against a physical fermion Green's function obtained 
via (^) using an identical number of sites in fig. ||b). 
The agreement, albeit qualitative, is rather satisfying. 
In Hb) the thin dashed line refers to the pseudo-fermion 
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FIG. 8. Comparison of spectra for a dimerized 20-site chain 
for various values of momentum fc = to tt and ti = t2 = Jo, 
5 = 0.048. (a) ED data, reproduced from ref.^'^ (b) solid 
line SCBA+RPA spectrum Ai{k,ij), dashed line SCBA-only 
spectrum j4JJ(fc,cij). 77 — 0.05. (a) and (b) are shifted relative 
to each other due to differing zeros of energy. 

spectrum. Obviously the similarity of this function to 
the ED result is less convincing thereby demonstrating 
the relevance of ground-state quantum fluctuations. The 
difference between the pseudo- and the physical-fermion 
spectrum is particularly evident in the incoherent part at 
wave vectors larger than tt/2. 

Closing this section we comment on the so-called 
dimerized i-model limit, i.e. ti = t2, Ji 3> J2 with 
Ji 0, and = J3 = 0. In this case the ground state of 
the spin system is a perfect dimer product-state while the 
only energy scale of the system is the hopping amplitude 
t = ti. The local density of states p{uj) for a single hole in 
one dimension is known to be p{lu) = l/\/uj'^ — 4^2-.which 
is independent of the particular spin-backgroundtH Cal- 
culating pdi^^Lj) — ^ y^t. ACT^°''(fc,a;) using the repre- 
sentation (^2[) and (^3[) we obtain the result shown in fig 
^. The SCBA-I-RPA approach shows some deviation from 
the exact result regarding the low-frequency regime and 
the overall band-width. The latter effect is reminiscentjjaf 
similar findings from analysis of the 2D AFM i-modelEj. 
The pure SCBA pseudo-fermion spectrum shows little 
resemblance with the exact result. 



VII. SPIN-LADDER SPECTRA 

In this section we turn to the case of J2 = i2 = which 
describes a two-leg spin ladder. In contrast to the defini- 
tion (|3^ ) of the fermions to account for the ID represen- 
tation, the unit cell remains identical to that of the lattice 
topology of fig. |l|. The pseudo-fermion operators on a 



7 



3, 
'a. 




FIG. 9. t- model DOS for iV = 256 and ti = t2. Solid, long 
dashed, and short dashed line refer to the SCBA+RPA, i.e. 
Pti(aj), SCBA-only, i.e. Pa(<^), and Brinkman-Rice result^^. 



rung can be labeled according to their (anti) bonding sym- 
metry and the reciprocal space is quasi two-dimensional 
with a wave vector k — (fc, ky) where ky = {n)0 which 
refers to the (anti)bonding state 



V2 ("^'^'-^ 



(34) 



Expressing G^(fc,a;) in terms of bk.k ,a renders the 
Green's function diagonal because the Hamiltonian con- 
serves parity and the operators (Q) exhibit a differ- 
ent signature with respect, to reflections at a plane per- 
pendicular to the rungsEil. The pseudo and physical 
fermion spectral functions are labeled accordingly, i.e. 
^o(7r) <j(^''^) ^^"^ ^o(ir) cr(^'^)' whcre the subscript 0(7r) 
refers to the (anti)bonding symmetry. 

Conservation of parity leads to an important con- 
straint regarding the pseudo-fermion self-energy. Since 
the spin-triplet has odd parity, scattering of a pseudo- 
fermion off a triplet can occur only via a transition be- 
tween the bonding and the antibonding fermionic states 
which have different parity. This implies a certain 'ro- 
bustness' of the bare pseudo-fermion bands of Hq of 
( pO| ) against renormalization by V because of the fi- 
nite energy gap 1t\ between the bare bands given by 
£7r(o)(fc) = +i3Cos(fc). Therefore, perturbative 

analysis of the single-hola-jdynamics on a laddec. using 
methods complementarytHlEj as well as relatecoa to 
that of this work have led to rather good agreement with 
numerical studies of finite system. 

In fig. |lo|a) the (anti)bonding spectra are shown for 
a case of strong intra-rung exchange Ji/ii = 5, Ja/ii = 

0. 5. The number of rungs is set to iV = 256 for which 
finite-size effects are negligible. The ground-state in this 
case is close to a pure singlet product-state with almost 
no quantum fiuctuations. The spin-gap is relatively large, 

1. e. A/ii sa 4.5. The spectrum clearly displays the rem- 
nants of the (anti)bonding bare bands and the intensity of 
the incoherent contributions to the spectrum are weak. 
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FIG. 10. Spectral function j4qj^j ^ of the bonding (0) and 
antibonding (tt) bands of the two-leg ladder for various values 
of momentum fc = to tt. = 256, r] = 0.1. (a) ti = t^, 
.h/ti = 5, JsAi = 0.5. (b) ti = t3, Ji/ii = 1, Ja/ti = 0.1. 



Next, fig. [lO|b), with Ji/ti = 1, Jz/ti = 0.1, refers a 
smaller gap which accounts for the extra incoherence of 
the spectrum and the less correspondence between the 
bare and the renormalized bands. Moreover the high en- 
ergy parts of the two bands intersect. 

Unfortunately, at the isotropic point, i.e. Ji = J3, the 
LHP approximation breaks down, with the spin-i. 
ing for J3 ^Ji/2. Yet, from series expansionc 
and DMRGES it is known, that A/ J w 0.5 in that case. 
To account for this deficiency of the LHP we keep the 
isotropic set of parameters in all of ( p0[]2l| ) however we 
adjust J3 in ([7[]Tl|) such as to result in A/Ji = 0.5. We 
believe that this modification provides an approximate 
description of the renormalization of Uk which occurs be- 
yond the LHP approach. 

Fig. |Tl|a) and b) display the case of isotropic exchange 
and hopping, i.e. = ta = Ji = J3. Here the renor- 
malization of the free particle bands is very strong. The 
bonding band at highest frequency, is very flat. In the an- 
tibonding band the intensity at highest energy is a com- 
posite excitation of states from the bonding band and a 
single triplet. Its relatively large intensity is a result of 
parity conservation, which inhibits intraband scattering. 
In addition a lower energy feature, reminiscent of the 
free antibonding band, can be observed at in the vicin- 
ity of fc = TT. Fig. [ll|b) shows a stretched frequency 
window with the two bands of lowest binding energy. In- 
terestingly, the maximum of the dispersion is found off 



8 




-5 5 1 

CO/tj CO/tj 



FIG. 11. Spectral function ^o{7r) o- of the bonding (0) 
and antibonding (tt) bands of the two-leg ladder for vari- 
ous values of momentum fc = to vr at the isotropic point 
ti — = Ji — Ji. N = 256, = 0.1. (a) complete spectrum 
(b) details at low binding-energy. 

the mcpienta fc = or tt. This is consistent vifith ED 
studiescj and series expansionsEEl. 

Analogous to the case of the Unear chain we contrast 
our diagrammatic calculation against ED analysis. In fig. 
|l2|a) a spectrjjxn reproduced from the work of Haas and 
collaboratorsEHI is compared to a physical-fermion spec- 
trum obtained from the SCBA-I-RPA approach. The lat- 
ter spectrum is shown in fig. [l2[b). It compares well with 
the ED result. Finally the dashed line in p^b) displays 
the SCBA, i.e. pseudo-fermion spectrum. The latter 
spectrum lacks the large intensity in the high energy re- 
gion of the antibonding band and the displacement of 
the maximum of the dispersion at low binding energy 
mentioned in the previous paragraph. This demonstrates 
that an inclusion of ground-state quantum-fluctuations is 
mandatory for a proper description of the single-hole dy- 
namics in this regime of parameters. 




-10.0 -6.0 -2.0 2.0 -6 -2 2 6 
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FIG. 12. Comparison of spectra for an A'^ = 8-rung two-leg 
spin-ladder for various values of momentum fc = to tt. 
ti — = Ji — J3, and 77 = 0.1. (a) ED data, reproduced from 
ref.^° (b) solid line SCBA-fRPA spectrum A%{k,uj), dashed 
line SCBA-only spectrum A%(k.,u}). For better visibility the 
SCBA-only spectra have been shifted by Au/t — 2 along the 
X, an arbitrary amount along the y axis. 

and the scattering of the holes by massive triplet-modes 
of the spin system. Due to the spin-gap we find, that 
the single-hole excitations are quasi-particle like though 
strongly renormalized. The impact of ground state quan- 
tum fluctuations of the dimer state on the hole-spectrum 
has been elucidated. Our results compares well with nu- 
merical diagonalization of finite systems, both for chains 
and ladders. We hope that our work will promote fu- 
ture experimental studies to obtain angular resolved pho- 
toemission spectra on the novel chain-and ladder-type 
transition-metal oxides in an energy window comparable 
to the magnetic exchange coupling. 
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VIII. CONCLUSION 

In conclusion we have investigated the single-hole ex- 
citations at half filling in dimerized and frustrated spin- 
chains and ladders in the spin-dimer phase. Based on 
an analytic approach we have provided a physically intu- 
itive picture of the single-hole dynamics as resulting from 
a combination of the dimerization of the spin background 
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